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Abstract 

A simple construction is presented, which generahses piecewise Hnear 
one-dimensional Markov maps to an arbitrary number of dimensions. The 
corresponding coupled map lattice, known as a simplicial mapping in the 
mathematical literature, allows for an analytical investigation. Especially 
the spin Hamiltonian, which is generated by the symbolic dynamics is 
accessible. As an example a formal relation between a globally coupled 
system and an Ising mean field model is established. The phase transition 
in the limit of infinite system size is analysed and analytical results are 
compared with numerical simulations. 
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1 Introduction 

The influence of chaotic motion on the dynamics of systems with a large number 
of relevant degrees of freedom is one of the central and unsolved issues of nonlinear 
dynamics. This problem has a strong impact on applications (cf. It is not 
very well understood from the analytical point of view, in despite of the enormous 
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amount of numerical results, which are available in the literature. This drawback 
comes to a certain extent from the lack of simple nontrivial models, which can 
be studied by elementary methods, and which allow for an investigation of such 
basic problems. The present publication intends to narrow this gap. 

Recalling the tremendous success, that low-dimensional time discrete dynam- 
ical systems had played in the understanding of low-dimensional chaos, systems 
of coupled maps have become a paradigm for the investigation of complex space 
time behaviour 0. They allow for both efficient numerical simulations and 
for analytical approaches. Especially it is possible to introduce a definite notion 
of space time chao^ 0, |^ . 

These analytical and even rigorous approaches are based on what is called the 
thermodynamic approach towards dynamical systems. Its main idea is based on 
the fact, that owing to a symbolic dynamics the stationary dynamical properties 
can be described in terms of the statistical mechanics of spin systems. Within 
such an approach one dimension of the spin lattice corresponds to the time axis 
of the dynamical system, whereas the other dimensions correspond to the spatial 
extension. Bifurcations in the dynamical system are refiected by phase transitions 
in the corresponding spin system. 

Such concepts are especially easy and elementary to apply within the context 
of one-dimensional piecewise linear Markov maps [§, 0, §|. Of course, since the 
dynamical system has no spatial extension the associated spin system is one- 
dimensional, and long range interactions are necessary to cause a phase transi- 
tion. Such interactions are typically related to a nontrivial grammar, singular 
expansion rates or the violation of transitivity in the dynamical system. 

Among spatially extended systems phase transitions may occur under much 
milder conditions since the associated spin lattice is at least two-dimensional. 
In fact, it is already a nontrivial task to establish a conventional high temper- 
ature phase, since the interactions which occur in the spin system are rather 
intricate. Nevertheless, it has been shown rigorously that certain weakly coupled 
maps admit correlations decaying exponentially in space and time, and that the 
corresponding spin system is in the paramagnetic phase. Hence it was suggested 
[0 that the occurrence of more complicated coherent patterns, usually observed 
in numerical simulations, are related to the phase transitions in the associated 
spin system. 

There are several examples available in the literature reporting phase transi- 
tions, and I will mention a few of them. For a certain class of coupled logistic 
mappings, where a space time mixing state in the weak coupling regime was 



proven to exist stable periodic patterns emerge beyond such a regime [p!T[] . 
Although these properties show rigorously that a phase transition occurs, noth- 
ing is known about its nature. Furthermore, the scenario resembles bifurcations 



""^As already pointed out by the authors, the terminology space time mixing is more 
appropriate. 
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in low-dimensional dynamical systems, since the spatial extension seems to play 
no crucial role for the transition. The same conclusion holds for the so called 



"peak crossing bifurcation" [12], a kind of boundary crisis which causes space 
time intermittency from a trivial non chaotic state. Finally the model proposed 
by Miller and Huse shows an Ising like transition, but only in two spatial 
dimensions. The model is quite different from that proposed here, since the local 
dynamics used by Miller and Huse mimics to a certain extent a single spin. In 
fact one-dimensional coupled map lattices of this type do not indicate a phase 
transition, but the real mechanism does not seem to be fully understood [|14| . 

Phase transitions associated with short range interactions in the symbolic dy- 
namics should display a finite size behaviour, since the transition is only possible 
in the thermodynamic limit, i.e. in the limit of large system size. Hence the cor- 
responding instability is expected to be related to a transient behaviour, where 
the length of the transients increase with the system size. Numerical simulations 
show such a behaviour (cf. |jT5|, |16[). Altogether these considerations indicate 
that it might be helpful to develop a model which can be handled analytically 
with moderate effort. 

I will propose a model which possesses a coupling quite different from the 
famous "diffusive" type coupling. It is not surprising that a model which can 
be handled analytically has a special type of spatial interaction too. For that 
reason section |^ is devoted to the introduction of the basic idea using two coupled 
maps only. We are lead to the discussion of simplicial mappings which are the 
natural generalisation of piecewise linear one-dimensional Markov maps. The 
construction will then be performed for an arbitrary number of maps in section 
1^. Although we will finally apply the results for a globally coupled system, the 
approach is certainly not restricted to that case as stated in the conclusion. 



2 Two— dimensional simplicial mapping 

Recalling that piecewise linear one-dimensional Markov maps can be handled 
analytically, it is tempting to construct a coupled map lattice with similar prop- 
erties. Particularly, I want to present a model which is by construction piecewise 
linear on its Markov partition, but allows for instabilities, that means phase tran- 
sitions, due to the spatial extension. However, before I dwell on this problem let 
me introduce in this section the basic ideas of this construction. For that purpose 
we restrict ourselves to the simple case of two coupled maps only. 

To begin with let us consider a system of two independent tent maps 

x^:l = {U x)t^ :=l-2|x(r)| . (1) 

^As already stressed in ^ this coupling does not mimic the real diffusion but is merely 
chosen as a simple and typical short range coupling. 
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Figure 1: Transformation properties of the uncoupled tent maps (|I]). A simplex 
and its image is indicated. 

Here the spatial index v takes the values or 1. Of course the Markov partition 
of the single map carries over as a direct product 

4(oM.):={(x(°),x«)|0<a('^)x(^)<l}, a^^) g {-1, 1} . (2) 

Each set (0) of the Markov partition is mapped to the full phase space by the 
dynamics (cf. fig.|I|). The grammar is trivial in the sense that all symbol sequences 
consisting of pairs (cr^^-*, cr*-^-*) are permitted by the dynamics. 

We are now going to introduce an interaction in such a way that the map 
remains piecewise linear on a suitable Markov partition and hence will allow 
for an analytical and elementary treatment. The strategy consists in fixing a 
suitable Markov partition, and then constructing the corresponding analytical 
form of the mapping. We start from the Markov partition of the uncoupled 
system and deform the partition appropriately (cf. fig.0). Consider a fixed set 
of the partition we intend to construct. On this set the two-dimensional map 
has to be affine. A general two-dimensional affine mapping has six adjustable 
parameters. Therefore by fixing three points in phase space and their images 
the affine mapping is determined uniquely. For that reason we cannot stick to 
the cubesQ, which constitute a Markov partition of the uncoupled system, but 
have to trace back to triangles, i.e. two-dimensional simplices. In fact, specifying 
a simplex and an image simplex fixes uniquely an affine transformation, which 
maps the former set onto the latter one. Therefore one introduces an interaction 
between the two tent maps, if one deforms the simplices in the Markov partition 

•^If one demands that the map is continuous and piecewise hnear on rectangles, then we are 
left with the uncoupled case. 
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Figure 2: Simplices of the Markov partition of two coupled maps obtained by a 
deformation of the uncoupled case. The shaded area indicates a simplex and its 
image. 

but requires that the images of these simplices coincide with the image simplices 
of the uncoupled case (cf. fig.@). However the deformation cannot be performed 
arbitrarily, since the image of each simplex has to consist of a union of simplices 
in order to guarantee the Markov property. Hence vertices on the boundary of 
phase space have to stay on the boundary, and the vertex on the diagonal has to 
stay on the diagonal. This simple geometrical prescription constitutes the basic 
idea for the construction of the coupled map lattice. In what follows a suitable 
notation will be introduced, in order to present the analytical expressions for the 
map and the associated partition function. The approach will be held general 
enough to capture the case of arbitrary dimension too. 

To label the simplices again we start from the uncoupled case. The simplices 
are obtained by dividing each square (0) into two simplices as indicated in fig.|l[ 
Each simplex has three vertices, the first being given by a "corner" of the phase 
space {a^'^\ o'^^^), the second being given by the intersection of the boundary and 
the coordinate axis, and the last one is located at the centre (0, 0). Corresponding 
to these vertices the following symbol is assigned to the two resulting simplices 

{(aW,a«),(0,a«),(0,0)} - 

{(aW,a«),(aW,0),(0,0)} ^ . (3) 

Alternatively this convention corresponds to the notation 

[(a(°U«);P] := {(x^, x^^)) | < a^^^^x^^ < < 1} , (4) 

if P denotes a permutation which is either Q, the transposition of the two 
numbers 0, 1 or J, the identity permutation. Within this notation the uncoupled 
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map To induces the following transitions, i.e. the grammar of the system 
[(^(0)^^(1))./] ^ [(l,l);g]u[(l,-l);g]U[(l,-l);/]U[(-l,-l);/] 
[(a(°),a«);Q] ^ [(1, 1); /] U [(-1, 1); Q] U [(-1, 1); /] U [(-1, -1); Q] .(5) 

As already mentioned the interaction between the tent maps is introduced, 
by shifting the vertices (±1, 0) and (0, ±1) along the boundary an amount, which 
we denote by t{±i,o) and t(o,±i) respectively. The shifted vertices read (±1, t(-i-i^o)) 
and (t(o,±i)) il)- In the same way the vertex on the diagonal (0, 0) is shifted to 
(^(0,0); ^(0,0))- The deformed simplex [{a^^\ a^-^^); P], to which we assign the same 
label, possesses the vertices 

The grammar and the condition that the map acts linearly on the simplices, 
define the dynamical system completely. Within this setting the numbers ^((t.o-') 
act as parameters for the map T. 

It is now possible to evaluate the partition function corresponding to the 
topological pressure. For that purpose we need the local expansion rate, i.e. the 
Jacobian. The latter is easily obtained from the ratio of the volume of a simplex 
and the volume of its image. Hence there is no need to write down the closed 
analytical formula for the coupled system at the moment. If one introduces the 
notation a = (cr'^°\ a^^^) and 

aW:=(0,a«), a« := ((t(°),0), := (0, 0) (7) 



then the volume of the simplex (0) is given by 
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Vol[K\aW),P] = la^(^)-t^,,,, a-^^^-Wu} 

= ^(l-^^(°)W.>)(i-^^«W.>) . (8) 

Since the volume of the image is always given by 2^/2 (cf. figj^) one obtains for 
the Jacobian the expression 

IdetDTl^,,, = 2^ (l - a^^%noyy' (l - a^^^^t^POu}) • (9) 

Performing the summation over all n-periodic symbol sequence^ [g_Q; Pq], . . . , 
[OLn-i'i Pn-i],[Qj3', Pq], ■ ■ ■ that are permitted by the grammar yields for the 
partition sum corresponding to the local expansion rate the expression 

^n(/3) = En2-^'(l-^?^\-.w) (l-^?^\-.ioa>) . (10) 



"'in fact, a slight ambiguity arises when introducing the partition function, since the corre- 
spondence between symbol sequences and orbits is not unique, if orbits on the boundaries of 
the partition are considered. 
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The summation in eq.(^) involves both a summation over symbol patterns {cr^*^ } 
and over sequences of permutations Pk which are related by the grammar. In 
order to read off the spin Hamiltonian from eq.(|10|) the summation with respect 
to the permutation sequences has to be performed. But as indicated in appendix 
1^ this sequence is almost uniquely determined by the symbol pattern, if surface 
terms are discarded. Hence expression ([T0|) may be understood as a sum over 
all periodic symbol patterns {(Tk"^}, with the sequence of permutations being 
determined by this pattern according to the grammar (^. Summarising, eq. ([TI1|) 
represents the partition sum of two coupled spin chains with the interaction being 
mediated by the deformation parameters t, i.e. by the coupling in the dynamical 
system. No phase transition occurs of course, since no long range interactions 
are involved. 



3 Ising type transition in a globally coupled model 

The considerations of the preceding section will be carried over to the case of 
L coupled maps, where L may be understood as the extension in the "spatial" 
direction of the dynamical system. 

We start again from L independent tent maps. Let 
specify the state of the system in the phase space [—1, 1]^. Each cube of the 
Markov partition, as usual labelled by G {—1, 1}^, is divided into L\ simplices 
according to (cf. eq.@) 

[ffl P] := {x I < < . . . < < 1} , (11) 

where P denotes a permutation of the L numbers 0, . . . , L — 1. The vertices of 
this simplex are given by 

. (12) 

Here a^i"' - ''^} denotes the phase space point which is obtained by replacing the 
entries at positions -P(O), . . . , P{v) in the symbol sequence g_ with zero. Uncoupled 
tent maps have the property, that the image TQ{g_^^^''"''^^) of a vertex is obtained, 
if each coordinate with value ±1 is replaced by —1, and each coordinate with 
value is replaced by +1. Using this rule the grammar corresponding to the 
partition (p!T| ) can be constructed, but we will not need these expressions in what 
follows. 

The interaction among the maps is introduced by deforming the simplices, 
i. e. by shifting the coordinates of the vertices (|12D which are zero by a constant 
amount. The shifted vertex is generated from a^'fO'---'''} by replacing every zero in 
the sequence by the same constant value t^p{o,...,u] , which of course is smaller than 
1 in modulus (cf. eqs.(^ and (^). The shift may differ from vertex to vertex 
as indicated by the subscript. For the shifted vertex the notation r^p{o,...,^} is 
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introduced. The construction of the coupled map lattice is completed by the 
condition that the image of the vertex is given by the image of the corresponding 
vertex of the uncoupled system T{t^p{o,...,v}) = To{a^^'^'■^^''^^). The last condition 
guarantees that the image of each simplex is given by a union of simplices like 
in the uncoupled case. Since a closed analytical formula of the coupled map 
lattice T will be of no use for us, I refrain from writing down such an expression. 
However for numerical purposes appendix |B| contains a simple recursive algorithm 
to evaluate T. 

In order to determine the Jacobian one needs the volume of the deformed 
simplices. The latter follows immediately from the vertices described above as 

Vol[a; P] = i- n (l - a^^%no,....y) . (13) 

Since the image of each simplex has the volume the Jacobian reads 

^-1 _i 
|detDT|[^.p] = 2^ n (l - ^''^'^^^{0,.....}) . (14) 

Hence we obtain for the partition sum 

/ n-l L-1 , . p 

UP) = E n 2"^^^ n (i - ^?^^^v.(o,....> ) , (15) 

k=0 1^=0 ^ -h J 

where again the summation has to be performed with respect to all n-periodic 
accessible symbol sequences [oio! -^0]; • • • ? \Q-n-i \ -Po]; [^o! -^0]; • • •• As already indi- 
cated in the preceding section the permutation sequence is (almost) uniquely 
determined by the two-dimensional symbol pattern {o'^^^^. Hence we may al- 
ternatively perform the summation with respect to all the spin variables, with 
the permutations being determined by the spin pattern. Then the correspond- 
ing Hamiltonian of the two-dimensional spin lattice can be read off from eq. ( [T5| ) 
directly. 

The structure of the Hamiltonian is however a little bit hidden since the 
permutations depend on the spin variables. To simplify the subsequent consid- 
erations we will refer in what follows to a special choice which corresponds to a 
globally coupled system, so that the permutations drop from the formulas. Let 
us choose the deformations in such a way, that on the one hand no interaction in 
the index fc, i.e. in the direction of time occurs, and on the other hand a global 
interaction between the spins in the index z/, i.e. an interaction in the spatial 
direction is realised. Then we are left with simple mean field Ising chains. 

n (1 - o^^^K,,,,...,.,)^ = exp (f3Hj:a^''^ + d ^ ^''"^cr^'^ + PEo] (16) 



8 



Confer appendix |C| for the corresponding explicit construction, which determines 
the deformation parameters t^p{o,...,,^} in terms of H and J and the normahsation 
constant £"0- Roughly speaking the parameter H, corresponding to the magnetic 
field in the spin system, describes the deviation of the map lattice from an inver- 
sion symmetric case, whereas the spin interaction J mediates the spatial coupling 
among the maps. Inspecting eqs.(^5|) and (^6|), the Hamiltonian of the full map 
lattice consists of non interacting spin chains. Since the interaction within a chain 
has infinite range a phase transition occurs at jSH = 0, (3 J = 1/2 if the thermody- 
namic limit L —>■ 00 has been performed. Hence our coupled map lattice displays 
a bifurcation, which is intimately related to the limit of large system size. Of 
course the model allows for an analytical calculation of those quantities, which 
are determined by the spin pattern, since the partition sum can be evaluated in 
the thermodynamic limit. 

Let us supplement the analytical results with numerical simulations of the 
full map lattice (cf. appendix and |C| for the numerical algorithm). Keep in 
mind however that the usual long time average, which will be denoted by (. . .) in 
the sequel corresponds to the choice = 1 in the context of dynamical systems. 
Obvious global quantities are given by the spatially averaged symbol sequence 
and its square 




The dependence of these averages on H and J is displayed in fig.^. The well 
known behaviour of the Ising mean field model is reproduced, but of course 
the results have been obtained here from the dynamics of the map lattice. One 
should however keep in mind that from the point of view of the dynamical system 
the quantities (|r^) look a little bit artificial and the actual interpretation of the 
dependence on the system parameters is difficult. Nevertheless the main feature, 
i. e. the second order phase transition clearly shows up at if = 0, J = 1/2, but 
finite size effects are superimposed. 

In order to demonstrate the influence of the system size on the transition more 
clearly, consider the zero field dependence of (S"^) on the interaction strength for 
different system size (fig.^). Since the phase transition sharpens if the thermo- 
dynamic limit is approached, it is obvious that the transition in the dynamical 
system is intimately related to the thermodynamic limit. Hence it is a character- 
istic feature of the high-dimensional dynamical system. 

Finally, let us have a look at the space time evolution in symbol space as well as 
in phase space. Consider first the evolution in an inversion symmetric situation, 
i.e. H = 0.0 (cf. fig. ^. Below the transition point, i. e. for weak coupling, a 
rather irregular state is observed. It corresponds to the space time chaotic regime 
mentioned in the introduction. Slightly above the transition point long range 
order develops, but a lot of "defects" are superimposed. The synchronisation 
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L=256 




Figure 3: Dependence of the spatially averaged symbol and its square on H and 
J for system size L = 256 and L = 128. The time average has been performed 
using a spatial uniformly distributed initial condition and 2000 iteration steps, 
after discarding a transient of length 200. 
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Figure 4: Dependence of (5^) on the interaction J for zero field H = and 
different system sizes (symbols). The solid line displays the analytical result in 
the thermodynamic limit L oo. 

becomes more perfect, if the coupling is further strengthened. Of course both 
phases are mixed, since no symmetry breaking field is present. With a symmetry 
breaking field the weak coupling regime behaves very similar to the zero field case 
(c.f. fig.^. But above the transition point one phase is dominant and a pure state 
develops if the coupling is increased. In summary, the phase transition shows up 
in the symbol representation as well as in the phase space coordinates, where 
apparently two different kinds of motion occur. 



4 Conclusion 

A class of piecewise linear coupled map lattices has been investigated which allows 
for an analytical approach. The key point of the construction consists in fixing a 
Markov partition in terms of simplices, imposing a grammar, and requiring the 
map to be linear on this partition. Such an approach allows for an elementary and 
analytical construction of the corresponding spin Hamiltonian. Phase transitions 
are possible at finite coupling in the limit of large system size, although the model 
is locally expanding and has the Markov property. 

The proposed construction may be understood as the natural extension of 
piecewise linear one-dimensional maps to more than one dimension. In fact, by 
a refinement of the partition a large class of functions, at least in finite dimen- 
sions, can be approximated, as stated by the simplicial approximation theorem. 
Furthermore, piecewise linear one-dimensional maps are able to approximate dy- 
namical features of low-dimensional chaos, as stated by the Ulam conjecture 



[0. Similar properties might hold also for simplicial mappings in connection 
with high-dimensional dynamics, but further clarification is certainly needed. 
The explicit evaluation of a phase transition has been performed for a global 
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Figure 5: Space time evolution of the symbols and of the phase space coor- 
dinate x^^'' for a system of size L = 256, field H = 0.0, and coupling strengths 
increasing from above to below J = 0.2,0.6, 1.0. Black/white pixels correspond 
to the values cr = —1/ + 1 and dark/light pixels to phase space coordinates close 
tox = -1/ + 1. 
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Figure 6: Space time evolution of the symbols and of the phase space coor- 
dinate x^^'' for a system of size L = 256, field H = 0.2, and coupling strengths 
increasing from above to below J = 0.2,0.6, 1.0. Black/white pixels correspond 
to the values cr = —1/ + 1 and dark/light pixels to phase space coordinates close 
tox = -1/ + 1. 
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coupling only. Although globally coupled systems are interesting in itself I have 
dealt with a globally coupled system for reasons of simplicity, since otherwise the 
exact computation of the phase transition may become much more involved. I 
believe that the principal aspects of this paper do not depend on this choice of 
coupling. Indeed an interaction in the time direction of the spin system may be 
introduced by tracing back to finer partitions, and a spatial short range coupling 
may be realised. Hence it seems possible to construct a model along these lines 
which is in some sense closer to a nearest neighbour coupled Ising model in two 
dimensions, i. e. a model with short range interactions in each direction but 
admitting a phase transition in the thermodynamic limit. Since the objective of 
the present publication was on the principal aspect I leave this problem for future 
work. 

In the treatment I have adopted a physicists approach, by considering sys- 
tems of arbitrary but finite size and performing the thermodynamic limit at the 
end of the calculation. Therefore it is not easy to decide rigorously whether 
the dynamical system, being defined by a geometrical construction, tends to a 
mathematically well defined limit. Nevertheless the discussion of the example 
supports such a conjecture. Then the phase transition like behaviour in the fi- 
nite model should be related to a transient dynamics whose duration increases 
with the system size. Such an analysis, which is of course desirable for a deeper 
understanding of the transition has not been accomplished yet. 

Models of the type investigated here, i. e. models with a large number of 
degrees of freedom showing a phase transition and allowing for an elementary 
analytical approach, may be helpful to understand several features of high- 
dimensional chaotic dynamics. Although the grammar, i. e. the structure of 
space time periodic orbits was fixed by hand, it is also possible to investigate the 
infiuence of pruning, by considering increasingly finer Markov partitions. These 
properties together with an analytical accessible system may have an impact on 
the understanding of expansions in terms of space time periodic orbits, which 
have recently attracted some interest. 
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A Two coupled maps 

In order to simplify our computations and not hide the main arguments behind 
extensive computations let us specialise to the case ^(0,0-) = t{a,o) — —c tanh J, 
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^(0,0) =0. I think it will become evident that this choice does not influence the 
generality of the arguments. 

The partition sum ([T0|) is taken with respect to all ra-periodic sequences 
{[njt'jPk])- Consider first those patterns which obey a^^^ ^ a^^^ for some index 
k E [0,n — 1]. Since the symbols {a^^\a^^^), that means the sets [a;/] U 
already determine a Markov partition, the pattern {cr'k^} specifies a phase space 
point, which in our case is a period-n orbit^ But this orbit does not hit a bound- 
ary of the simplices and especially not the diagonal of the phase space, since at 
least one pair of symbols has different entries. The validity of this statement is 
obvious for the case of noninteracting maps. It carries over to the general case, 
since the deformation of the simplices does not destroy topological properties in 
the phase space. Hence the sequence of permutations Pjf^ is determined uniquely 
by the periodic orbit, that means by the pattern of symbols {<y^k^}. 

Consider instead a'^^ = cr^^"* for all k, i.e. a synchronised orbit. Then the 
corresponding phase space point is located on the diagonal, i. e. on the boundary 
of the simplices. From the grammar (j^) we have that the predecessor of a 
permutation Pk+i changes according to the symbol cr^k+i^ ^^"^ obeys 

Pfc = Pfc+iif4+i'^ = -l Pfc = gPfc+i if = 1 . (18) 

The condition that the sequence Pk is n-periodic, Pq = Pn, requires that the finite 
sequence {a^^^\ . . . , (T^ll^) contains an even number of symbols +1. Furthermore 
there exist two corresponding permutation sequences according to eq.(|l^) and the 
choice Pq = I,Q. 

Then the partition sum (|I0|) splits into two parts, one part containing the 
summation over non synchronised symbol patterns, and the second part contain- 
ing the summation over synchronised symbol patterns with an even number of 
+ 1 symbols and two complementary permutation sequences. If we supplement 
the first part with the missing synchronised symbol sequences]^ then eq. ([T0| ) reads 

UP) = E n2-^4i-a?(V(o>)7i-^?^Vio.>'' 



k ' k 



„{0)_„{i) k=0 ^ -k / \ -k J 



(T, ' =(t; ' ,even 

k k 



E n2-^^(l-o-?^V(o>) (l-af^«V.ou>) (19) 



^ For the particular case cr^'^'' = 1, tr^^^ = — 1 such an argument does not apply. But this 
single point does not matter in the thermodynamic limit. 

^In the following an arbitrary but henceforth fixed construction rule for the permutation 
sequences is imposed. 
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where the index even/odd means, that the summation is performed with respect 
to sequences containing an even/odd number of +1 symbols. Now in each sum the 
sequence of permutations is uniquely determined by the symbol pattern {cl^^}- 
By our choice of deformations the relation 

1 - ^^''\.^oyy (l - -?^^^t^^.ou>)' = (1 + tanh J)^ (20) 

holds, and it is exactly this property which simplifies the subsequent consid- 
erations. Furthermore we perform in the second and third sum of eq.(19) an 
unconstrained summation over ak, < k < n — 2, choosing an-i in such a way 
that the constraint of an even or odd number of +1 values is guaranteed. Then 
we get 

n— 1 „ 

UP) = E n2-^'(l + 4°^4'^tanhj)' (21) 

n-2 

+ n2"'^(l + tanhJ)^{(l + tanhJ)''-(l + tanhJ)''} . 

ao,...,a„-2 k=0 

The second contribution vanishes identicall}^ Now the Hamiltonian can be read 
off directly from the partition sum 

n-1 

-/57^ = (^-Jai°Vi^^ + lncoshJ + ln4) . (22) 

k=0 

We end up with a simple pair interaction in the "spatial" direction and no inter- 
action in the "time" direction. 



B L— dimensional simplicial mapping 

For a given set of deformation parameters t^p{o,. ..,,,} a recursive algorithm will be 
developed, which allows for the calculation of the image point T{x). Let us 
assume that x is contained in the interior of a simplex, which from the point of 
view of a numerical evaluation causes no restriction for the continuous mapping T. 
Since the simplices are deformed, it is a nontrivial task to determine the simplex, 
that contains a given phase space point x. However using the fact that simplices 
are convex, a sequence of projections is adopted to determine the corresponding 
set, the simplicial coordinates and hence the image point (cf. fig.|3). 

Let X[o] = X denote the phase space point, and Tq := t^p{o,...,l-i} = (to, to, ... , to) 
with to ■= t^p{o,...,L-i} = t(o,...,o) the central vertex in the cube [—1, l]'^. By our 

^In general, that means for maps which are not symmetric with respect to the origin, this 
contribution seems to be a surface term, i.e. it can be neglected compared to the first sum if 
the exponential n-dependence is considered. 
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Figure 7: Diagrammatic view of the recursive projections onto lower-dimensional 
faces of the simplex, which are used to determine the simplex containing the phase 
space point x. The box indicates a quadrant [0, 1]"^, denote the vertices of the 
simplex, and xj^j the sequence of the projected points. 



construction this vertex is common to all the simplices and hence is a vertex of 
the simplex which contains the phase space point x. Consider the semi-infinite 

I[o]- 

at a point denoted by x^ij 



ray emanating from Tq through xroj. It intersects the surface of the cube [—1,1]-'' 



^[1] = soX[o] + (1 - so)to . (23) 
Here sq denotes the number which makes the expression 

-So := — Tin 7r-\ *■ (24) 

|„(fco) ^ I (ko) ^ ^ ^ 

F[o] ''Ol ^[0] ''O 

minimal among all components k. The component ko for which the minimum is 
attained determines the surface of the cube where the intersection occurs by the 
condition x^'^°^ — a^^°\ with a^^°^ being given by 

= sign(a:[J]°) - to) ■ (25) 

That surface contains one face of the simplex we are looking for. If we define 

P(L — 1) := fco) then the central vertex within this surface reads T_i :— T(^p{o,...,l-2}, 
and it leads to the second vertex of the simplex we are looking for. Now one 
repeats the same reasoning using the points X]^^ and Ti, but with the dimension 
of the problem being reduced by one. 
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Hence one continues by induction. Determine fc^i ^ {P{L — 1), ... , P{L — /i)} 
among those components, which have not been fixed yet, such that 

■= \J>^^)_f I " ^^^^ 
is minimaL Here := t^p{o,...,L-i-M}. Then 

a^'^^ ■.= signix[';f -t,) (27) 

and P{L — 1 — fi) := kfj_ determine the next components of the symbol sequence 
and the permutation respectively. 

= + (1 - s^)t^ (28) 

yields the new intersection point, using the abbreviation = t^p{o,...,l-i-ij.}. This 
iteration terminates exactly after L steps, since we end up with a corner point of 
the cube [—1,1]^. It yields the simplex [a; P] which contains the point x. 

In addition we can read off the simplicial coordinates < < 1 which are 
defined by 

x=T.>^^r,, E^^ = l ■ (29) 
^l=o ^i=o 

These coordinates are particularly useful for the evaluation of the image point 
T{x). Recalling that the map T is linear on the simplex, and that the images of 
the vertices are obtained, by replacing the components with values ±1 by — 1 
and the components with values by +1, we get 

(T(x))^(^) = t>^.{nr,)f'\ E K{nr^)Y^'^ = ±X,- E A, 

/i=0 fi=u+l fi=0 ii=i/+l 

V 

= -1 + 2EA^ . (30) 

Furthermore we can express the simplicial coordinates in terms of the quantities 
(P^D, which have already been determined in the course of the recursion. As we 
prove at the end of this appendix 

t^, = ^-if[^)i (31) 

ti=0 \ii=0 J 

holds for < z/ < L — 1. Hence the simple recursion scheme described above to- 
gether with eqs.(^) and (|31| ) complete the algorithm for the numerical evaluation 
of the map lattice. 
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We are left with verifying eg. (pT]) . Again an induction will be used. First the 
component P{L — 1) of eg. (|29|) yields 

x^(^-i) = Vo + (1 - AoV^^"-^) . (32) 
Then using the definition p5| ) one obtains 

1 1 to 



(33) 



which by virtue of eg.(piD leads to the result (^TJ) for z/ = 0. Suppose that (PI|) 
holds for V < p. Evaluation of eg. (|29|) for the component P{L — p — 2) yields 

^p{L-p-2) ^ ^^^^ ^ ^ _ ^ Ap+itp+i + (1 - Ao - . . . - Ap+Oa^^^-''-^) . (34) 
On the other hand, we obtain from the recursion relation (L 



PiL-p-2) _ , . ( PiL-p-2) (1 - So)to (1 - Si)ti 

(1 - Sp)tp tp+i \ 



Sq . . . Sp Sq . . . Sp 



(35) 



If one now inserts eg.(0) into eg.(p5D and repeatedly uses eg. (|3l|) for < p one 

gets 

P(L-p-2) 

= (SO...S,) f(l-Ao-...-A,+i)a^(^-^-^) + A,+iVi ^ 



\ Sq . . . Sp J 

= (i_Ao-...-A,+i)(.o...3p)(a^(^-''-^)-Vi) . (36) 
But then from the definition (^) we have 

- - - = (so . . . Sp) I p(L-p-2) ^ : ~ p(L~p-t) . I (37) 

and the result follows from the definition (^). 

C Spatial mean field coupling 

Using the abbreviation 

Su := E ^^"^ (38) 
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a sequence of functions gu{Su) will be constructed which obeys 



L-l 



2J 

T 



-S, + h) +J + Eo = ^ln(l-a('^)tanh(7,(5,)) (39) 

with a constant Eq being independent of the spin variables. Since the left hand 
side is invariant with respect to a permutation of the spin numbers the result 
(|TBp follows then immediately if one defines 



t^p{o,...,^} := tanhgfj. 



E 



cr 



(40) 



^M^P{0,..,i/} 



Applying the trivial identity 

/(a) = a[/(+l) - /(-l)]/2 + [/(+1) + /(-l)]/2 (41) 

one obtains 

In (l - (t(^) tanh (7.(5,)) = -a('')(7,(S,) - KiS,) (42) 

if the abbreviation ht,{x) := In cosh5f,^(x) is introduced. Observing 5*1, = a^'^'^^^ + 
Si,+i and applying eq.(^) again to the second term on the right hand side one 
gets 

K{S.) (43) 



a 



2 2 



If we repeat these steps for the last term on the right hand side and recall that 
by definition 5^-1 = 0, one ends up with 



L-2-V 

E - 



1 n+l 
(n+i/+l) ^ 



+ 



n=0 



)L-l-^ 



2n+l 



E 

A;=0 



fe=0 

L - 1 - 
k 



n 



h^{n + 1 - 2k + Sn+i+u] 



K{L-l-v-2k) 



(44) 



Insert eqs. (^^ and (^) into eq. (pQ]) and equate the coefficients of a^'^) on both 
sides to obtain finally 



2 7 y \ p 



k=0 



p-l 
k 



p-l 
k-l 



K-p{P -2k + Sy 



(45) 

This system of equations defines recursively the functions gy{SS) for increasing 
index v. Since the independent variable is discrete and takes only a finite number 
of values, the numerical evaluation is quite straightforward. It should be remarked 
that a simple asymptotic evaluation for large system size L is possible too. Such 
an evaluation indicates that gu{Su) remains uniformly bounded in the limit of 
infinite system size. 
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